Theory of radiative and rare B decays 

Gino Isidori 

INFN, Laboratori Nazionali di Frascati, Via E. Fermi 40, 1-00044 Frascati, Italy 

Abstract. We present a concise theoretical overview of radiative and rare B decays mediated by flavour-changing neutral- 
current transitions of the type b — » s(d)y and b — > s(d)M. 



INTRODUCTION 

Thanks to the efforts of B factories, the exploration of 
the mechanism of quark-flavour mixing is now enter- 
ing a new interesting era. The precise measurements 
of mixing-induced CP violation and tree-level allowed 
semileptonic transition have provided an important con- 
sistency check of the SM, and a precise determination 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The 
next goal is to understand if there is still room for new 
physics (NP) or, more precisely, if there is still room 
for new sources of flavour symmetry breaking close to 
the electroweak scale. From this perspective, radiative 
and rare B decays mediated by flavour-changing neutral 
current (FCNC) amplitudes represent a fundamental tool 
(see e.g. Ref. (J). 

Beside the experimental sensitivity, the conditions 
which allow to perform significant NP searches in rare 
decays can be summarized as follows: i) decay amplitude 
dominated by electroweak dynamics, and thus enhanced 
sensitivity to non-standard contributions; ii) small theo- 
retical error within the SM, or good control of both per- 
turbative and non-perturbative corrections. In the rest of 
this talk we shall analyze at which level these conditions 
are satisfied in various decay modes. 



INCLUSIVE FCNC B DECAYS 

Inclusive rare B decays such as B — > X s y, B — > X s t + t~ 
and B — » X s vv are the natural framework for high- 
precision studies of FCNCs in the AB = 1 sector |2J]. Per- 
turbative QCD and heavy-quark expansion form a solid 
theoretical framework to describe these processes: inclu- 
sive hadronic rates are related to those of free b quarks, 
calculable in perturbation theory, by means of a system- 
atic expansion in inverse powers of the fe-quark mass. 

The starting point of the perturbative partonic calcula- 
tion is the determination of a low-energy effective Hamil- 



tonian, renormalized at a scale fi = G(m^), obtained by 
integrating out the heavy degrees of freedom of the the- 
ory. For b —> s transitions -within the SM- this can be 
written as 

J%g = --l=V*V tb £ Q(M)Gi + h.c. (1) 
V2 i=1 

where (2t...6 are four-quark operators, Q% is the chromo- 
magnetic operator and 
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are the leading FCNC electroweak operators. Within 
the SM, the coefficients of all the operators in Eq. @ 
receive a large non-decoupling contribution from top- 
quark loops at the electroweak scale. But the m t depen- 
dence is not the same for the four operators, reflecting 
a different SU (2)i-breaking structure, which can be af- 
fected in a rather different way by new-physics contribu- 
tions nn. 

The calculation of the rare decay rates then involves 
three distinct steps: i) the determination of the initial 
conditions of the Wilson coefficients at the electroweak 
scale; ii) the evolution by means of renormalization- 
group equations (RGEs) of the C\ down to [i = G(m\,)\ 
iii) the evaluation of the hadronic matrix elements of the 
effective operators at jj. = ffimy), including both pertur- 
bative and non-perturbative QCD corrections. Each of 
the three steps must be taken to matching orders of ac- 
curacy in powers of the strong coupling constants a s and 
of the large logs generated by the RGE running. The in- 
teresting short-distance (electroweak) dynamics that we 





FIGURE 1. Representative diagrams for the mixing of four- 
quark operators into Q-j (left) and Qg (right). 



would like to test enters only in the first step; the follow- 
ing two steps are fundamental ingredients to reduce and 
control the theoretical error. 

The first two steps (initial conditions and RGEs) are 
process independent and are common also to exclusive 
modes. Nonetheless, the organization of the leading-log 
(LL) series is not the same for the three underlying 
partonic processes, or the four operators in Eq. 10: 

b — > sy. Here only Qi has a non-vanishing matrix el- 
ement at the tree level. The large logarithms generated 
by mixing of four-quark operators into Qi (see Fig. ^ 
play a very important role and enhance the partonic rate 
by a factor of almost three 01 . Since this mixing van- 
ishes at the one-loop level, a full treatment of QCD cor- 
rections beyond lowest order is a rather non-trivial task. 
This has been achieved already a few years ago, thanks 
to the joint effort of many authors (see e.g. Ref. ylEj] and 
references therein), and is nowadays a rather mature sub- 
ject. All the ingredients of the partonic calculation have 
been cross-checked by more than one group. In partic- 
ular, very recently an independent confirmation of the 
three-loop mixing of Qi and 2i...6 1 7] - till few months 
ago the only piece of the calculation performed by one 
group only - has been presented @] . 

b — > s£ + £~ . The three operators with non-vanishing 
tree-level matrix elements are Qi, Qg and £>io- Similarly 
to Qi, QCD corrections are very important also for Qg. 
Since Qg mixes with four-quark operators already at the 
one-loop level (see Fig. 0, the organization of the LL 
series for b — > s£ + £~ is different than in b — > sy: the NLL 
level is much simpler (no three-loop mixing involved), 
but less precise. An accuracy below the 10% level on 
the decay rate (a precision similar to the NLL level in 
b — > sy) is reached here only with a NNLL calculation. 
All the missing ingredients to reach this goal has finally 
become available. In particular, NNLL initial conditions 
and anomalous dimension matrix can be found in and 
respectively. 

It is worth to stress that the impact of QCD correc- 
tions is very limited in the axial-current operator 2io- 
This operator does not mix with four-quark operators 
and is completely dominated by short-distance contri- 



butions. Together with Q v , (2io belongs to the theoret- 
ically clean G(Gp) hard-GIM-protected part of the ef- 
fective Hamiltonian. Thus observables more sensitive to 
<2iO, such as the forward-backward (FB) lepton asymme- 
try in b — > s£ + £~ , have a reduced QCD uncertainty and a 
stronger sensitivity to possible non-standard phenomena. 

b — > svv . In this case only Q v is involved. Similarly 
t° Gio> QCD corrections play a very minor role since 
there is no mixing with four-quark operators. As a result, 
the only non-trivial step of the perturbative calculation 
for b — > svv decays is the determination of the initial 
condition of C v at the electroweak scale: this is known 
with a precision around 1% within the SM I Toll . 

These two processes-independent steps of the calcula- 
tion, namely the determination of the effective Hamilto- 
nian renormalized at a low scale fi = 0(m^\ can easily 
be transferred from the b — > s case to the b — > d one. The 
only difference is the richer CKM structure of the b — > d 
Hamiltonian, with two independent non-negligible terms 

(v; d v tb mdv: d v uh ). 

The situation is very different for the last step of the 
calculation, namely the evaluation of the hadronic ma- 
trix elements. The latter strongly depend on the specific 
process and the specific observable we are interested in 
(e.g. fully inclusive rate or differential distribution). In 
the following we shall review the results of this step (and 
thus the final numerical predictions) for some of the most 
interesting b — > s observables. 



B-*Xsy [the most effective "NP killer"] 

The inclusive B — > X s y rate is the most precise and 
clean short-distance information that we have, at present, 
on AB = 1 FCNCs. Combining the precise measurements 
by ALEPH, BaBar, Belle and CLEO, the world average 
reads jH 



SS{B^X s yf^ = (3.34 ±0.38) x 10" 



(3) 



On the theory side, the NLL partonic calculation of 
the matrix elements, performed first in Ref. lUl for 
the leading terms, has recently been cross-checked and 
completed in Ref. lll3ll . Perturbative corrections due to 
higher-order electroweak effects have also been analyzed 
(see Ref. 1 14] and references therein). 

Non-perturbative 1 /ntf, corrections are well under con- 
trol in the total rate. In particular, ^(l/m&) corrections 
vanish in the ratio T(B -> X s y)/T(B -> X c £v), and the 
0{l/m 2 b ) ones are known and amount to few per cent 
11511 . Also non-perturbative effects associated to charm- 
quark loop s have been estimated and found to be very 
small Il6l fl7l Hill . The most serious problem of non- 
perturbative origin is related to the (unavoidable) experi- 
mental cut in the photon energy spectrum that prevents 



the measurement from being fully inclusive [l^, fl^l - 
With the present cut by CLEO E y > 2.0 GeV EU, this 
uncertainty is smaller but non-negligible with respect to 
the error of the perturbative calculation. The latter is 
around 10% and its main source is the uncertainty in the 
ratio m c /mj, that enters through charm-quark loops Bill . 

According to the detailed analysis of theoretical errors 
presented in Ref. Ell , the SM expectation is 



(4) 



B — ► X s £ + £ [the present frontier] 



?(B^X s y) SM = (3.73 ±0.30) x 10~ 4 



in good agreement with Eq. Q. It must be stressed 
that the overall scale dependence is very small: for jx £ 
[nti,/2,2mi,] the central value moves by about 1%. The 
error in Eq. (0} is an educated guess about the size of 
possible NNLL terms. In particular, the largest source of 
uncertainty is obtained by the variation of ra c (jU)/»jP ole 
for ji £ [m c ,tn.j,]. A critical discussion about the error in 
Eq. 10}, with alternative more conservative estimates, can 
be found in Ref. 1 22] . 

The comparison between theory and experiments in 
S3(B — > X s y) is a great success of the SM and has led 
us to derive many significant bounds on possible new- 
physics scenarios. For instance, the 3§(B — > X s y) con- 
straint is one of the main obstacles to build consis- 
tent models that predict a sizable difference between 
J^ce{B — * QKs) and s^ce{B — * yK$). This constraint can 
be avoided (see e.g. Ref. 1 23 ] and references there in), 
but the resulting models require a considerable amount 
of fine tuning. By far more natural are the so-called MFV 
models 0), where £?cp(B — > <j)K s ) w £#cp{B — * yKs) 
and deviations from the SM in 38(B — > X s y) do not ex- 
ceed the 10%-30% level |4]. Improved measurements of 
S$(B — ► X s y) are certainly useful to further constrain this 
possibility. However, since the experimental error has 
reached the level of the theoretical one, it will be very 
difficult to clearly identify possible deviations from the 
SM, if any, in this observable. 

Hopes to detect new-physics signals are still open 
through the CP-violating asymmetry 



T{B^X s y)~T{B^X s y) 
T{B^X s y)+T(B^X s y) 



(5) 



Ar CP (B^X s y) 



This is expected to be below 1% within the SM 1,221124 
but could easily reach values beyond the SM, 

even in the absence of large effects in the total B — ► X s y 
rate. This is indeed one of the main expectations in 
models with sizable NP effects in £/qp(B — >• <pK$). The 
present measurement of AFcp(fi — > X s y) is consistent 
with zero 11 ill , but the sensitivity is still one order of 
magnitude above the SM level. 



Both Belle I25H and BaBar 1 26] have recently an- 
nounced a clear evidence (« 5cr) of the B — > X s £ + £~ de- 
cay. The two results are compatible and are both based 
on a semi-inclusive analysis (the hadronic system is re- 
constructed from a kaon plus to 4 pions, with at most 
one 7t°). Their combination I Till 

33{B — + X s £ + £~) exp = (6.2± l.l+};f) x 10~ 6 . (6) 

represents a very useful new piece of information about 
AB = 1 FCNCs, with considerable margin of improve- 
ment in the near future. 

In principle, these decays offer a phenomenology 
readier than B — > X s y, with more than one interesting ob- 
servable. The joint effort of several groups has recently 
allowed to evaluate at the NNLL level all the matrix el- 
ement necessary for the two main kinematical distribu- 
tions: the dilepton spectrum 1 2*71 EH l29ll and the lepton 
FB asymmetry I30ll3l"ll . 

In addition to the non-perturbative corrections due to 
the finite b quark mass, B — ► X s £ + £" transitions suffer 
from specific non-perturbative effects due to long-lived 
cc intermediate states (B — > X s cc — > X s £ + £~). The heavy- 
quark expansion, which allow to evaluate the Aqcd 
terms, is rapidly convergent and leads to small correc- 
tions for sufficiently inclusive observables ll32l l33ll . A 
consistent treatment of the second type of effects re- 
quires kinematical cuts in order to avoid the large non- 
perturbative background of the narrow cc resonances 
(see Fig. 0. These two requirements are somehow in 
conflict 1281 13311 : nonetheless, we can identify two per- 
turbative windows, defined by: 



1 



M: 



e+t- 



S [1 GeV 2 , 6 GeV 2 ] (low) , 
(high) , 



q 2 > 14.4 GeV 2 

where reliable predictions can be performed [28}. It is 
worth to emphasize that the two regions have comple- 
mentary virtues and disadvantages: 

• Virtues of the low-q 2 region: reliable q 2 spectrum; 
small I /nib corrections; sensitivity to the interfer- 
ence of C7 and C9; high rate. 

■ Disadvantages of the low-q 2 region: difficult to per- 
form a fully inclusive measurement (severe cuts on 
the dilepton energy and/or the hadronic invariant 
mass); long-distance effects due to processes of the 
type B -> yx, -> X s +X'£+£- not fully under con- 
trol; non-negligible scale and m c dependence. 

• Virtues of the high-q 2 region: negligible scale and 
m c dependence due to the strong sensitivity to 
|Cio| 2 ; easier to perform a fully inclusive mea- 
surement (small hadronic invariant mass); negligi- 



ble long-distance effects of the type B 

X s +X'£+£-. 




FIGURE 2. Dilepton spectrum of the inclusive B — * X s e + e 
decay within the SM. Vertical axis: d.3S(B — > X s e + e~)/dq 2 in 
units of 10~ 7 x GeV~ 2 ; horizontal axis: q 2 in GeV 2 .The dotted 
line denotes the NNLL pure perturbative result, the full line 
includes an estimates of the non-perturbative cc effects 1 28 ] . 



We stress that this prediction is already saturated by ir- 
reducible theoretical errors and, contrary to the results in 
Q, is very difficult to improve it further. 

As anticipated, some of the most interesting short- 
distance tests in S$(B — > X s £ + £~) decays can be per- 
formed by means of the FB asymmetry of the dilepton 
distribution: 



1 



d@(B^X s t+t-)/dq 2 J -i 
d 2 38(B -> X s £+£-) 



d cos 9, 



dq 2 d cos 0( 



sgn(cos d e ) , (8) 



where 0( is the angle between l + and B momenta in the 
dilepton centre-of-mass frame. Here the SM predict a 
zero for s = ql/m 2 = 0. 162 ± 0.008 HHHl: a very pre- 
cise prediction which could easily be modified beyond 
the SM, even in absence of significant non-standard ef- 
fects on the total rate. 



• Disadvantages of the high-q 2 region: q 2 spectrum 
not reliable (only the integrated rate can be pre- 
dicted); sizeable 1 /m/, corrections (effective expan- 
sion in 1 / (nib — \/<7min) GUl); low rate. 

Given this situation, we believe that future experiments 
should try to measure the branching ratios in both regions 
and report separately the two results. These two mea- 
surements are indeed affected by different systematic un- 
certainties (of theoretical nature) and provide a different 
short-distance information. The NNLL SM predictions 
for the two clean windows I2S1, 



(B^X s £ + r)Z 

M 

gh 



i(B^x s e+r)^ 



= (1.63 ± 0.20) x 10 
= (4.04 ± 0.78) x 10 



-6 



-7 



(7) 



are still affected by a considerable error; however, in both 
cases the uncertainty is mainly of parametric nature and 
could be substantially improved in the future. In particu- 
lar, the large error in the high-g 2 region is mainly due to 
the uncertainty in the relation between the physical q 2 in- 
terval and the corresponding interval for the partonic cal- 
culation (i.e. the uncertainty in the relation between nij, 
and the physical hadron mass), which can be improved 
with better data on charged-current semileptonic modes. 
According to the recent analysis of Ref. 1 29], both results 
in Q should be decreased by sw 4% to take into account 
the leading electroweak corrections. 

The two results in Q cannot be directly confronted 
with 0, which includes an extrapolation to the full q 2 
spectrum. The updated SM expectation for this extrapo- 
lated branching ratio is (4.6±0.8) x 10~ 6 1 28] (in agree- 
ment with the previous estimate of Ref. |34]), and it is 
consistent with the present experimental world average. 



EXCLUSIVE MODES 

On general grounds, theoretical predictions for exclu- 
sive FCNC decays are more difficult. The simplest cases 
are processes with at most one hadron in the final state. 
Here there has been a substantial progress in the last few 
years, both by means of analytic approaches 1 35] and by 
means of Lattice QCD 1 36], but still the overall theo- 
retical uncertainty is around 20% at the amplitude level. 
The largest source of uncertainty is typically the normal- 
ization of the hadronic form factors, an error that can 
be substantially reduced in appropriate ratios or differ- 
ential distributions. These type of observables become 
particularly interesting in channels where, because of ir- 
reducible experimental problems, the short-distance am- 
plitude cannot be extracted from corresponding inclusive 
modes. Two of such examples are the ratio 



R 7 (P/K*) = 



,<%{B^py) 
38{B -> K*y) 



(9) 



and the normalized FB asymmetry in B — > K*l + £ . 

The ratio Ry(p /K*) is one of the most promising tool 
to extract short-distance properties about the b — > sy 
amplitude. On the experimental side, the combination 
of the bounds on charged and neutral channels, in the 
isospin limit, leads to R r (p/K*) < 0.047 at 90% C.L. 
13711 - On the theory side, the B — > Vy amplitudes which 
determine this ratio have been analyzed beyond naive 
factorization by several authors 138.13^. 1401. Within the 
SM one can write 



RJp/K* 



Vtd 
V, s 



( m l-m 2 rho ) 



{m\-m\,y 



-C 2 (l-Afl), (10) 



Wpy/k'Y) / 



/ R(PY/KY) 
New lattice result 



FIGURE 3. Bounds on the p-T] plane from the BaBar bound 
onRy(p/K*): the two curves correspond to different values of 
? 139]. 



where £ denotes the ratio of the form factors at = 
in the mt, — > 00 limit, and A/? the additional SU(3) (and 
isospin) breaking due to and 0(oc s ) corrections. 

The largest source of uncertainty is £: according to the 
light-cone sum rule estimate £ = 0.76 ±0.10 139I1 . the 
present experimental bound on Ry(p /K*) is about twice 
the SM expectation. However, preliminary Lattice results 
indicate a larger value, £ = 0.91 ±0.08 |41], which 
would imply more stringent bounds on | V tc t /V ts | . The two 
curves in Fig.[5]summarizes the present status. 

While the inclusive B — > X s £ + £~ rate is already ac- 
cessible at the B factories, a differential study of the 
inclusive FB asymmetry is beyond their present and 
near-future reach. More accessible from the experimen- 
tal point of view is the FB asymmetry in B — > K*£ + £~ 
(defined as in Eq. (JSJl with X s — > K*). Assuming that the 
leptonic current has only a vector (V) or axial-vector (A) 
structure (as in the SM), the FB asymmetry provides a 
direct measure of the A—V interference. Indeed, at the 
lowest-order one can write 



^ FB (<7 2 )-Re^C 



'10 



1 n sS 
— 2 C 9 - 



2 ra^Cv 

r{q 

m B 



where r{q 2 ) is an appropriate ratio of B — > K* vector and 
tensor form factors [42]. There are three main features 
of this observable that provide a clear and independent 
short-distance information: 1) The position of the zero 
of .(^fb(<? 2 ) m the low-g 2 region (see Fig. @J 1 42] . As 
shown by means of a full NLO calculation |40], the ex- 
perimental measurement of q\ could allow a determina- 
tion of C7/C9 at the 10% level. 2) The sign of £?m(q 2 ) 
around the zero. This is fixed unambiguously in terms 
of the relative sign of C\q and Cg |43]: within the SM 
one expects s^^{q 2 > <?o) > for \B) = \bd) mesons. 
3) The relation £?[B] m (q 2 ) = -£/[B] FB (q 2 ). This fol- 
lows from the CP-odd structure of jz^b an d holds at the 
10 -3 level within the SM |43], where C10 has a negligible 
CP-violating phase. 
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FIGURE 4. Zero of the forward-backward asymmetry in 
B~ -* K*-£+£- at LO and NLO. The band reflects all theo- 
retical uncertainties from parameters and scale dependence 
combined 1 4fll . 



B s ,d — * [the future frontier] 

The purely leptonic decays constitute a special case 
among exclusive transitions. Within the SM only the 
axial-current operator, Q\o, induces a non-vanishing con- 
tribution to these decays. As a result, the short-distance 
contribution is not diluted by the mixing with four-quark 
operators. Moreover, the hadronic matrix element in- 
volved is the simplest we can consider, namely the B- 
meson decay constant 



(0\qYnY5b\B q (p)) = i Pll f Bq 



(11) 



Reliable estimates of fs d and fs s are obtained at present 
from lattice calculations and in the future it will be pos- 
sible to cross-check these results by means of the B + — > 
£ + V rate. Modulo the determination of fs q , the theoreti- 
cal cleanliness of B s j — > £ + £~ decays is comparable to 
that of the golden modes K L — * n°vv and B — * X s d vv. 

The price to pay for this theoretically-clean amplitude 
is a strong helicity suppression for £ = ji (and £ = e), 
or the channels with the best experimental signature. 
Employing the full NLO expression of C10 itioll . we can 
write 



(B, 



. M + M -)SM = 3.1x10-9 



\Vrs\ 



fa 



0.21 GeV 

S8(B S 



1.6 ps 

SM 



KB, 



■ n+n- 



0.04 

m,{rn t ) 
166 GeV 

: 215 . 



3.12 



The corresponding B c ( modes are both suppressed by an 
additional factor \V td /V ts \ 2 = (4.0 ±0.8) x 10~ 2 . The 
present experimental bound closest to SM expectations 
is the one obtained by CDF on B s — > jJ. + fi~ fTllH : 



38{B S -> < 9.5 x 10" 



(95% CL) 



which is still very far from the SM level. The latter will 
certainly not be reached before the LHC era. 

As emphasized in the recent literature G^lEfiH^ . the 
purely leptonic decays of B s and B c / mesons are excel- 
lent probes of several new-physics models and, partic- 
ularly, of scalar FCNCs. Scalar FCNC operators, such 
as 1>rSlP.rHl, are present within the SM but are ab- 
solutely negligible because of the smallness of down- 
type Yukawa couplings. On the other hand, these am- 
plitudes could be non-negligible in models with an ex- 
tended Higgs sector. In particular, within the MSSM, 
where two Higgs doublets are coupled separately to up- 
and down-type quarks, a strong enhancement of scalar 
FCNCs can occur at large tan/3 = v u /v ( i |45|. This effect 
is very small in non-helicity-suppressed B decays and 
in K decays (because of the small Yukawa couplings), 
but could enhance B — ► £ + £~ rates by orders of mag- 
nitude. As pointed out in Ref. 14811 . 0(100) enhance- 
ments in 33{B — ► £ + £~) correspond to breaking 
of universality in 58(B — > Kii + n~) vs. 88 {B — > Ke + e~). 
Therefore, the present search for B — > £ + £~ at CDF is 
already quite interesting, even if the sensitivity is well 
above the SM level. In a long-term perspective, the dis- 
covery of such processes is definitely one of the most 
interesting items in the B-physics program of hadron col- 
liders. 



CONCLUSIONS 

Rare FCNC decays of B mesons provide a unique op- 
portunity to perform high-precision studies of quark- 
flavour mixing. The B — > X s y rate, where both exper- 
imental and theoretical errors have reached a compa- 
rable level around 10%, represents the highest peak in 
our present knowledge of FCNCs. The lack of devia- 
tions from SM expectations in Y(B — > X s y) should not 
discourage the measurement of other clean and indepen- 
dent FCNC observables, such as the forward-backward 
asymmetry in B — > X s £ + t~ or the B — > £ + £~ rates. Even 
if new physics will first be discovered elsewhere, the ex- 
perimental study of these theoretically-clean observables 
would still be very useful to investigate the flavour struc- 
ture of any new-physics scenario. 
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